Comments on nonlinear viscosity and Grad's moment method 
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It is shown that the steady unidirectional flow with vanishing heat flux considered by B. C. Eu [Phys. Rev. 
E 65, 031202 (2002)], and earlier by Uribe and Garcia-Colfn [Phys. Rev. E 60, 4052 (1999)], is inconsistent 
with the exact conservation laws of mass, momentum, and energy. The inconsistency does not lie in the assumed 
symmetry properties of the flow but in the stationarity assumption. The unsteady problem is considered and its 
solution from the Boltzmann equation for Maxwell molecules is given. 
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In a recent paper 1 1], Eu analyzed a steady unidirectional 
flow at uniform temperature and derived the equations for the 
stress tensor elements from the Boltzmann equation by means 
of Grad's moment method. These are essentially the same 
state and the same method as those considered by Uribe and 
Garcia-Colm 1 2], except that the transversal velocity gradients 
were assumed to vanish in Ref. |2], while they are included in 
the analysis of Ref. 1 1 ]. The major aim of this Comment is to 
show that the steady unidirectional flow at uniform tempera- 
ture of Refs. (J 01 is inconsistent with the exact macroscopic 
conservation equations. 

For a dilute monatomic gas, the macroscopic balance equa- 
tions expressing the conservation of mass, momentum, and 
energy are |3l||4|] 
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where Dt = (3( + u • V is the material time derivative, n is the 
local number density, u is the local flow velocity, T is the local 
temperature, m is the mass of a particle, kB is the Boltzmann 
constant, q is the heat flux vector, and P is the pressure (or 
stress) tensor The flow considered in Refs. 1 1, 2] is character- 
ized by the following properties (not necessarily independent): 
(a) it is a unidirectional flow, i.e., u(r) = Ux{y)% where x is 
the unit vector along the flow direction; (b) the temperature 
is uniform, VT — 0; (c) the heat flux vanishes, q 0; (d) 
the pressure tensor is uniform; and (e) the state is stationary, 
i.e., dt 0. Let me consider first the geometrical properties 
(a)-(d) separately from the stationarity assumption (e). Appli- 
cation of assumptions (a)-(d) on the exact balance equations 
0-0 yields 
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dtn + (nu^) = 0, 
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Are Eqs. @-(|6j consistent with a steady state? Equation Q 
shows that the flow velocity is stationary if and only if the flow 
is incompressible, i.e., if dux/dx = 0. In that particular case, 
Eq. (0} is consistent with a stationary density if and only if 
the density is uniform as well. But, even if dux/dx = and 
n — const, energy balance equation (|6jl implies that the tem- 
perature cannot be stationary but monotonically increases with 
time due to viscous heating effects (note that Pxydux/ dy < 
and Pxzdux/dz < because of physical reasons). There- 
fore, the steady state assumption (e) is incompatible with as- 
sumptions (a)-(d), except in the trivial case Vux — 0, i.e., at 
equilibrium. The unsteady Boltzmann equation for the incom- 
pressible unidirectional flow with dux/dx — dux/dz — 0, 
dux/dy — const, usually referred to as uniform shear flow or 
homoenergetic simple shear flow, has been solved exactly for 
arbitrary shear rates in the case of Maxwell molecules |5, 6, 7]. 
An analogous solution has been obtained in the case of the BGK 
model kinetic equation for more general interactions 1 8 , 9 , 10]. 

In the case of a compressible flow in the absence of transver- 
sal gradients, i.e., Vuj; ||x, Eq. Q shows that the flow velocity 
is necessarily unsteady. According to continuity equation @, 
it is still mathematically possible that dtn = if the product 
nux is uniform, i.e., n{x)ux{x,t) = K{t). Insertion of this 
condition into Eq.Q yields i^^^i^(i) n^^n'{x) — —A^^, 
where the dot denotes a time derivative, the prime denotes 
a spatial derivative, and A is a constant. The solution to 
these equations is simply n{x) = A/{x — xq), Ux{x,t) = 
a{x — xo)/ (1 + at), where Xq and a are constants. This math- 
ematical solution is unphysical unless the problem is restricted 
to the half domain x > xq (x < xq) if A > (A < 0). But 
even in that case the existence of a nonuniform density is in 
conflict with the uniformity assumptions (b) and (d) because 
in a dilute gas n = p/ksT, where p = ifr P is the hydrostatic 
pressure. In summary, assumptions (a)-(d) do not contradict 
the conservation laws Q-© in the compressible flow with 
dux/dx 7^ if and only if the three hydrodynamic quanti- 
ties (density, flow velocity, and temperature) are unsteady, so 
assumption (e) is again incompatible with (a)-(d). 
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Strictly speaking, assumptions (b) and (d), i.e., uniform 
temperature and pressure tensor, were not explicitly stated in 
Ref. 101 ■ Actually, assumption (b) is implicit in (c); otherwise, 
one would have a thermal gradient that does not produce any 
heat flux, what is at odds with the second principle of thermo- 
dynamics. As for assumption (d), it was replaced in Ref. 0] 
by a weaker one: (d') the irreversible part of the pressure ten- 
sor is uniform, namely d{Pxx — p)/dx ~ 0. From a physical 
point of view it seems difficult to reconcile a uniform normal 
stress difference P^x — P with a nonuniform hydrostatic pres- 
sure p. In any case, let me drop conditions (b) and (d) for the 
moment and prove that conditions (a), (c), (d' ), and (e) are also 
inconsistent with the conservation laws, except at equilibrium. 
Application of (a), (c), and (e) on Eqs. ([3-<|3} gives 
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mnux 



dx 



Ux + -^Pa 
ax 
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where, as in Ref. |2], the case Vu^^Hx has been considered. 
Equations can be easily integrated to get 



nUx = riouo = const, 



= Pq + mrtgUQ = const, 



(10) 
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where the subscript denotes quantities evaluated at some 
reference point x = xq. So far, assumption (d') has not beeen 
used. This condition implies that Pyy — p = Pq — po = const, 
so that, according to Eq. (^J, 



p + mnu^ — Pa + ^7^f^o^^o ~ const. 
Insertion of Eqs. il 1> and ( I13t into Eq. il2\ yields 



(13) 



3 /I 

2^0 + + mnouo ( -uq - 2ux 



{ux - uo) = 0. (14) 



Both solutions of this quadratic equation are constants (but 
the physical one is Ux = uq). This closes the proof that as- 
sumptions (a), (c), (d'), and (e) are not consistent with the 
conservation equations except in the trivial case Ux = const. 

On the other hand, there is nothing wrong with assumptions 
(a)-(d) in the unsteady case. What is then the right form of n 
and Ux if Vm^ ||x? Since assumptions (b) and (d) imply that n 
is uniform, Eq. states that —dux/dx = n~^h = —K{t), 
so Ux{x,t) = K{t){x — Xq). Substitution into Eq. (|5|i gives 



K~^K = —1, so we finally have n{t) = no/(l + at), 
Ux{x, t) — a{x — xo) / {1 + at). This simple flow is known as 
homoenergetic extension |0, fllll . Again, the unsteady Boltz- 
mann equation can be solved exactly for arbitrary values of the 
constant control parameter a in the case of Maxwell molecules 
ll^QllQ^I, as well as with Grad's method lITsIl and in the case 
of the BGK model kinetic equation fl?) for more general in- 
teractions. A situation where a transversal velocity gradient 
duy/dx coexists with a longitudinal one dux/dx has been 
studied by Galkin flTll . 

It is illustrative to recall the application of Grad's method to 
the unsteady unidirectional flow at uniform temperature with 
no transversal gradients 1 13|. We have seen above that conser- 
vation of mass and momentum imply that = no/(l + at) 
and (t) = dux{x,t)/dx — an{t)/nQ. Positive values of 
the longitudinal deformation rate ^x represent some sort of 
"explosion" (or expansion) flow, while negative values repre- 
sent an "implosion" (or condensation) flow The energy 
balance equation (|6} becomes 



dtp 



:Px 



(15) 



In Grad's method, this equation is coupled to the (approximate) 
evolution equation for the stress element Pxx u3' 



dt {Pxx -P) + yPxx ~P + '^Pxx j 7x = (Pxx - P) , 

(16) 

where nonlinear terms have been neglected on the right-hand 
side and /i = p/r/j^s is an effective collision frequency, ryNs be- 
ing the Navier-Stokes shear viscosity. Equation il6\ is equiv- 
alent to Eq. (21) of Ref. |2], except that the time derivative 
is absent in the latter Without the time derivative operator, 
however, Eq. il6\ cannot be made consistent with Eq. ( I15t . 
After a transient regime, the system reaches a generalized 
hydrodynamic regime with Pij{t) = p{t)P*j{'Jx{t))^ where 
7*(t) = jx{t)/fj,{t) is the reduced longitudinal deformation 
rate. In general, from Eqs. il5i and il6i one gets a nonlinear 
first-order ordinary differential equation for P*^ (7*) |12, 1^. 
In the special case of Maxwell molecules (i.e fj, (x n), 7* (but 
not jx) is independent of time so one gets from Eqs. il5\ and 
( I16t an algebraic quadratic equation whose physical solution 

is p;, = 3 - 2p;^ with 



P* h*) 



8l*x 




:l*x + ^rf - 1 + 27^ 



(17) 



It turns out that this result for Maxwell molecules goes beyond 
the scope of Grad's method since it can be exactly derived from 
the Boltzmann equation 16, 11. 121 . 

To the best of my knowledge, Grad's method has not been 
applied yet to the unsteady unidirectional flow at uniform tem- 
perature with both longitudinal and transversal gradients. For 
this flow, an analysis of Eqs. (0} and (|5|l shows that, in addi- 
tion to n{t) = no/(l + at) and ^x = dux/dx = an{t)/nQ, 
one generally has shear rates jyx = dux/dy — a2n{t)/no 
and -fzx = dux/dz = a^n{t)/nQ, where a, 02, and 03 are in- 
dependent constants. Standard application of Grad's method 



yields the following evolution equation for the elements of the 
stress tensor: 



dtPr, 



dux 
' dx 



3 

-E 

k=l 



(18) 

Taking the trace in this equation one recovers the energy con- 
servation equation (|6j. It is worth noting that Eq. ( I18> can 
again be obtained exactly from the Boltzmann equation in the 
case of Maxwell molecules. For this interaction potential, the 
reduced rates 7* = -/^/fJ., 7^^ = lyx/f^, and 7*^ = j^^/fJ- 
are constants and Eq. ( I18> yields a coupled set of algebraic 



equations for P*j 
can choose 7*^. = 
that case, Eq. ilHi yields 



Pij/p- Without loss of generality 
7* , so that P* = P*, and P* = 
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FIG. 1: Plot of the normal pressure element Pyy relative to the hy- 
drostatic pressure p, as derived in Refs. Lln|2j|, versus the longitudinal 
deformation rate in the range —0.4 < 7J < for (a) zero shear rate 
(lyx = 0) and (b) 7*^ — i^The thin solid lines correspond to the 
linear approximation, Eq. <21> . while the dashed line in case (a) cor- 
responds to the nonlinear approximation, Eq. 122\ . The thick solid 
lines represent the exact results for Maxwell molecules. 



These two equations include as particular cases the homoen- 
ergetic extension flow (7*^ — 0), in which case Eq. ( I17t is 
recovered, as well as the uniform shear flow (7* = 0), where 

4: yy ' 



xy — 4^ yy \^ ^yy) ^yy solution of the cubic 

equation 47;^2^*,3 ^3 fl 



p* 



yy " \" " yy)' 
The inconsistency of assumption (e) on the stationarity of 

the flow geometrically characterized by (a)-(d) [or (a), (c), 

and (d')] manifests itself in the results obtained in Refs. |1] 

and 1121 for the stress tensor elements from Grad's method. 

Assuming cylindrical symmetry (Pxy — Pxz, Pyy = Pzz), Eu 

gets the following expression for the normal pressure element 

Rjy (neglecting nonlinear terms in the collisional integrals) 



II: 



P, 



yy 



, (l + 27*)(l + 37*) 
'(l + 27*)(3 + 77*)+47* 



(21) 



According to Eq. J21> . Pyy becomes negative at least in the in- 
terval — I < 7^ < — |, regardless of the value of 7*^. This is 
unphysical because the diagonal elements of the pressure ten- 
sor are positive definite quantities. Setting 7*^ — 0, Eq. i21\ 
reduces to the result derived by Uribe and Garci'a-Colm in 
the linear approximation |j2|]. When the nonlinear terms are 
included, they get |^ 
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P* =8 + — 7* 
yy ' o IX 



94 49 
1 H H 1* 



(22) 



The 



According to this expression, Pyy < if 7* < — ^• 
prediction of negative values of Pyy can be observed in Fig.Q] 
which shows the ratio Pyy/p given by Eqs. i21\ and i22\ in 
the range —0.4 < 7* < for vanishing shear rate (7*^ — 0) 
and for 7*^ — ^ . Comparison with the exact results illi for 
-f*x = and Eqs. (O and (|20} for j*^ in the case of 
Maxwell molecules shows that the predictions i2l\ and i22\ 



are only valid in the Navier-Stokes domain of small gradients, 
where Pj,j,/p « I + I7*. It is worth noting that the application 
to the unidirectional flow of a rheological theory by Eu 1 1 5 1 also 
yields unphysical negative values for the diagonal elements of 
the pressure tensor 1 16]. 

In Ref. 1 1], Eu claims that Grad's moment method is not 
thermodynamically consistent 1 17|. Actually, Grad's method 
is not but an approximate scheme for (partially) solving the 
hierarchy of moment equations stemming from the Boltzmann 
equation. The point I want to emphasize is that the physical 
inconsistency of the equations for the stress tensor elements 
derived in Refs. 1 1, 2] does not lie in the use of Grad's method 
(even if nonlinear terms are neglected), but in the wrong ansatz 
about the stationarity of the flow. In fact, as said before, the ap- 
plication of Grad's method to the unidirectional flow at uniform 
temperature for Maxwell molecules gives the same evolution 
equations for the stress elements as the Boltzmann equation. 
Therefore, at least in this instance, Grad's method is free from 
any thermodynamic inconsistency. 

It might be argued that assumptions (a)-(e) [or (a), (c), (d'), 
and (e)] are used in Refs. 1 1, 2] only as a tool to derive rheo- 
logical constitutive equations relating the irreversible part of 
the stress tensor to the velocity gradients in a nonlinear way 
by means of Grad's method. Such constitutive equations could 
then be applied to the conservation equations Q-0 regard- 
less of whether the flow is steady or not, whether the pressure 
is uniform or not, etc. However, it is doubtful that a constitu- 
tive equation derived from assumptions incompatible with the 
conservation laws can be acceptable beyond the Navier-Stokes 
regime, as Fig.^illustrates. 

Before closing this paper, let me comment on a couple of 
remarks made in Ref. 1 1 ] which are not directly related to the 
discussion made so far First, Eu states that "the shear viscosity 
is impossible to define" in the absence of transversal velocity 
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gradients and, consequently, "it is impossible to measure a 
shear viscosity without shearing the fluid." This is mislead- 
ing. To clarify this point, take the Navier-Stokes constitutive 
equations, namely ||41, 

P„=,,„_^,(gl + ||-|v.u*«). (23) 

In the special case of a unidirectional flow u = Ur^x, Eq. j23> 
yields 

Pxy ^ -mS-Q;^ , ^yPxx - Pyy) = -mS.-Q^- (24) 

Thus, to Navier-Stokes order, the response of the shear stress 
Pxy to a shear rate du^/dy is the same as the response of 
the normal stress difference {Pxx — Pyy)/'^ to a longitudinal 
deformation rate dux/dx. As a consequence, the Navier- 
Stokes shear viscosity can be measured from the normal stress 
difference, even in the absence of shearing (dux/dy = 0). 

The second point refers to Eu's claim 1 1 1 that the velocity dis- 
tribution function /(r, v; t) obeying the Boltzmann equation 
must always depend on the three spatial coordinates despite the 
fact that the hydrodynamic variables may depend on one space 
coordinate only, e.g., n{x, t), u{x, t), T{x, t). In support of 
this claim, Eu recalls that "even if the fluid particle moves one 
dimensionally in its hydrodynamic configuration space, it does 
not mean that the molecules making up the fluid particle and 
contained in the elementary volume of the hydrodynamic con- 
figuration space [. . . ] should be moving one dimensionally." 
While the quoted sentence is entirely correct, Eu's conclusion, 
namely that one cannot have /(r, v, t) — f{x, v, t), does not 
apply to the Boltzmann velocity distribution function / (r, v, t) 
but to the microscopic one-body distribution function defined 



by 

N 

F{r, v,t)=Y, - r,{t))S{v - v,(i)), (25) 

i=l 

where {ri(t),i = 1, . . . , TV} and {vi(t), i = 1, . . . ,7V} are the 
sets of positions and velocities of the particles of the system at 
time t. Actually, the velocity distribution function /(r, v, t) is 
the average of F{r, v, <), 

/(r, V, t) = (F(r, V, t)) ^ J dT i^(r, v, t)p(r), (26) 

where p(r) is the probability density or ensemble for the initial 
state and the integration is carried out over all the points F of the 
phase space. While in a given microscopic realization of the 
system i^(r, v, t) is a highly nonuniform function, its statistical 
average /(r, v, i) has a much smoother spatial dependence. 
In particular, it can depend on one coordinate only or it can 
even be uniform (e.g., at equilibrium). Of course, the fact 
that the hydrodynamic fields have a one-dimensional spatial 
dependence does not necessarily mean that the same holds to /, 
but there is nothing wrong if one restricts oneself to solutions to 
the Boltzmann equation with the same symmetry properties as 
the hydrodynamic fields. In fact, the so-called normal solutions 
are those that depend on space and time through afunctional 
dependence on the hydrodynamic fields |4J. 
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